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fh'r  reps rt  represents  «n  atteant  to  bring  together  from  the  fields 
*:  -  .  iynn.-.  I  cs ,  aerodynamics,  a  no  mathea-atics  oil  of  the  elements  '  i-  the  theory 

•rr.i'-r.ytr.;'  Method  of  Charac  ter  inties  in  establishing  the  fundamental  relatione 
in  the  uuf-rcenic  flew  of  a  perfect  gas.  References  1,  2,  ana  $  have  been  used 
ter.;-.  ively  in  gathering  -this  material. 

A  brief  description  of  eziroir.g  characteristics  flow  programs  is  also 
included.  More  detailed  information  concerning  specific  programs  is  available  in 
'-he  "arena rut  Data- Processing  Department. - 

II.  IMTROIUCITON 


These  programs  have  been  used  in  the  design  of  high  Mach  number, 
Internal  -  external,  variable  geometry  inlets  and  of  high  expansion  ratio,  vari¬ 
able  geometry  plug  nozzles.  Test  .’•'••suits  from  both  inlets  end  nozzles  have  shown 
,c  *i  >  r.  1 2  *  -j  y  o  ^  n  £  v  L  t»ri  p  re  diet  v  d  r.  03  c  re  « I  c.  a  1  pc  r  Co  r  in  2  nc  e  • 


Extensive  use  has  also  been  mace  of  these  programs  in  analytical  in¬ 
vestigations  of  flow .phenomena  in  supersonic  combust  ion  processes. 

Presently  existing  programs  will  handle  isentropic,  isent’nalpic, 
Irriaticnal  flew  '  for  both  two  -d  i  mens  ional  and  exi symmetric  designs.  Calorific' 
gas  imperfections  are  accounted  for  to  the  extent  indicated  in  the  body  of  this 
report.  Plow  conditions  behind  a  conical  shock"  can  also  be  computed. 

'Mork'is  proceeding  on  other  programs  which  consider  the  effect  of 
such  flew  phenomena  as  chemical  reactions  in  the  flow,  rotational  flow,  shock 
reflections  from  a  boundary,  boununry  lay  r,  and,  shock  wave-boundary  layer  inter¬ 
actions'. 

III.  .SASIC  THERMO  S'/IIAKICS 

'  _  A.  General  Relations  •!  -  - .  - 


In  order  tc  'present  the  analysis  for  supersonic  flow,  of  a  perfect 
t  is  first  necessary  to  write  down  a  relation  among  the  state  variables  P, 
(f  of  perfect  gases.  Such  a  relation  is  called  on  equation  of  state. 


;h“  ecuation  of  state  of  a  nerfvet  gas  is 


v  -  <7  np 
*  *  .  ' 
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An  alternate  formulation  Is 

PV  -  RT 


(3) 


Where 

V  *  Volume  per  unit  mass 


A  distinction  should  be  made  betveen  thermally  perfect  and  calorific- 
ally  perfect  gas-  A  thermally  perfect  gas  Is  one  which  obeys  Equations  (l)  or 
(2).  A  calorlfically  perfect  gas  la  one  whose  specific  heats  are  constant.  The 
geo  assumed  for  the  existing  characteristic  programs  is  thermally  perfect  and 
calorlfically  imperfect.  Calorific  imperfections  are  accounted  for  by  providing 
a  curve  of  the  ratio  of  specific  heats  (  P)  **  a  function  of  temperature. 

'The  perfect  gas  law  may  be  derived  from  the  following  :echanical 

models 

1.  The  gaa  is  assumed  to  ccnsist  of  a  set  of  point  masses  enclosed 
in  a  cubical  box. 

2.  The  point  masses  are  in  random  motion,  colliding  with  each 
other  and  tne  walls  of  the  box  in  perfectly  eleotic  fashion. 

From  this  we  can  conclude  th*t  the  change  of  momentum  per 
unit  time  at  each  wall  of  the  box  is  the  same. 

3.  Define  temperature  T  such  that 

T  —  2  Vr2,  or  | XT  -  |  i*Vr2, 

Where 

Vr  -  Root  mean  square  velocity  of  the  point  masses 

m  <•  Mass  of  each  point  mass 
k  ■  Boltzmann  constant 

This  establishes  a  relation  between  (molecular)  kinetic  energy  and 
temperature.  Real  gases  can  of  course  not  be  exactly  represented  by  this  model. 
However,  the  behavior  of  real  gases  can  be  adequately  described  by  Equations  (l) 
or  (2)  for  host  applications.  In  many  problems,  the  assumption  of  perfect  gas 
behavior  is  far  less  restrictive  than  certain  other  assumptions  which  are  made. 

B.  Equlpsrtltlon  of  Energy 

In  Item  3  above,  temperature  T  is  related  to  the  kinetic  energy  of 
translation  of  point  masses.  If  we  no v  consider  molecules  rather  than  point 
masses,  we  see  that  there  are  other  modes  of  motion  other  than  translational. 
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cf  inert  gooes,  which  contribute  to  the  total  energy  of  the 
can  be  related  to  temperature.  These  modes  of  motion  are 
jfid  vibrational  motion. 


of  a  diatomic  molecule,  the  total  energy 


oUcr.»ir.g: 


Translational  energy,  whereby  the  center  cf  mass  cf  the 
tr.oi ecu] p  ,r.ov»s  n  the  x,  Y,  and  Z  directions,  {three 
"3©gr£«£S  oT  freedom) 

RctatSonal  one* ivy,  vh*?reby  th*?  molecule  rotates  about  the 
center  cf  r.as,; ,  as  a  rigid  rotator,  { tve  degrees  cf 
freedom) 

Vibrational  energy,  whereby  the  atoms  of  the  molecule  vi¬ 
brate  along  ‘-he  axis,  or  line  connecting  the  atoms,  (one 
cc-gree  of  freedom) .  Ac  a  first  approximation,  we  con 
assume'  that  simnle  ha rr.cn ic  motion  is  executed. 


..n 


general,  the  degrees 
Translation: 

Pot  Fit  5  cr.  i  >-■  fer 


rf  freedom  are  allocated  as  follows: 


a  linear  molecule,  y  for  a  nonlinear 


molecule,  0  fer  a  monatomic  molecule, 

Vibration:  an  -  w  for  a  linear  molecule,  in  -  6  for  a 

nonlinear  molecule,  0  fer  a  ncrarttcmic  molecule. 


Humber  of  atoms  in  the  molecule 

The  Principle  of  Zaui partition  cf  Knergy  asserts  that  with  each 
-T  freedom  is  associated  energy  1/2  KT. 

In  dealing  with  the  vibrational  energy,  wo  assume  that  the  potential 
:  on  the  eve rage  equals  the  Kinetic  energy.  This  also  follows  from  the 
.tier*  cf  s iicole  harmonic  action. 


at  Constant  Volume 


"r.m.  the  above,  ve  car.  deduce  un  expression  fer  the  internal,  or 
energy  tf  a  gas.  If  there  ere  P  molecules  cf  gas  per  unit  macs,  then 


-f  f4.  .  -}•  K-  .  ,  j 
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Where 

E  *  Internal  energy  of  the  gas  per  unit  mass 

N  «  Number  of  degrees  of  freedom,  plus  the  vibrational  degrees  of  freedom 
for  potential  energy. 

(Actually,  l ho  Internal  energy  of  a  real  gas  also  depends  very  slightly  on  the 
volume  and  temperature.)  (We  define  a  calorlflcally  perfect  gas  os  one  for  which 

*  -  V(t)  ). 

Oy 

g 

Define  Cv  •  ~  «  Specific  heat  at  conatar.t  v'-lurne  lor  a  perfect  gas. 

Thua,  for  a  diatomic  molecule,  the  number  of  degreea  of  freedom  is 
6,  and  H  Is  6  +  1  *  7,  so  E  *  -2-  RT  and  Cv  «  R. 

D.  First  Law  of  Thermodynamics 


The  flrat  lav  of  thermodynamlca  states  that,  for  a  closed  system 
containing  a  gas  and  a  quasistatic  proceaa. 


dE  -  Sq^  ♦  £w. 


(«♦) 


Where 

£  Qq  *  Quantity  of  heat  exchanged  between  the  system  and  the  surroundings 
£>W  -  Amount  of  work  done  on  the  gas  by  the  surroundings 
dE  ■  Change  in  internal  energy  which  results  from  £q  and  S  W. 

If  the  internal  Energy  of  the  gas  is  a  function  of  the  state  vari¬ 
ables  only,  dE  la  a  perfect  differential,  but  $W  and  §Qq  are  not  necessarily 
perfect.  The  symbol  $  la  used  to  Indicate  that  these  quantities  are  not  func¬ 
tions  of  only  the  initial  and  final  states,  but  also  depend  upon  Intermediate 
conditions. 

We  now  obtain  a  relation  between  Cv  and  C_,  the  specific  heat  at 
constant  pressure.  The  Internal  energy  of  a  gas  depends  upon  P,  T,  and  P  ,  or 
P,  T,  and  V.  From  the  equation  of  state,  which  we  will  take  to  be  the  perfect 
gas  law,  we  can  express  any  one  of  P,  T,  or  V  in  terms  of  the  other  two.  Thus 
E  is  an  explicit  function  of  any  two  of  P,  V,  or  T. 
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r'r'i::  the  first  law, 


$  v  +  Sq0  =  "  (•g-m)  d?  +  (--gy)  dv 

V  T  r 

If  there  id  no  volume  change,  tiae®  and  (— ffi-)  «• 

=  dQ0- 


-  01 T,  wc  obtain 


(  ~  *c  \  -  e  «-( -•?  -  ?  -') 

•  Lv  'TTjT*  • 

J  ‘  V  5  T  V 

Rewriting  the  first  law,  ond  using  the  perfect  gae  result  that 


d(pV) 

AT 


CY  **  >  *  CV 


'of ine 


.  -  *vC  . 


E.  Er 


The  enthalpy,  or  total  heat  of  a  gaa  la  defined  to  be 


h  «.  PV  ♦  E 

■;sin,  this  is  referred  to  unit  mass.  From  the  results  for  a  perfect  gaa. 


h  -  rt  ♦  c^t_  .v:  , 


Here  tho  reference  level  Is  taXcn  ter®#;  .  h  *>  0  wjfeeo  T  ■  0. 
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We  must  now  distinguish  between  reversible  and  irreversible  proc- 
esses.  If  when  the  gas  is  performing  work,  or  is  being  worked  upon,  conditions 
exist  such  that  the  gas  is  not  in  equilibrium  with  the  surroundings,  the  gas  will 
undergo  certain  net  internal  accelerations.  In  order  to  return  to  equilibrium, 
this  type  of  motion  must  be  dissipated  as  heat  in  the  gas.  This  amount  of  wasted 
heat  is  added  to  the  term  SQq  in  t*»  eiotencnt  of  the  first  law  for  reversible 
proceecr- ,  Is  give  the  new  term  ^Q: 


$>%  * 


Define  the  entropy  (S)  by 


Ss- 


Note  that  SO^iss  is  not  which  is  exchanged  between  the  system 

and  the  surroundings,  but  heat  which  is  generated  in  the  system  itself  by  dissi¬ 
pative  processes.  In  the  case  of  the  system  doing  work  adlabatically  (  bQo  “  0), 
the  dissipation  due  to  friction  serves  to  reduce  the  amount  of  work  capable  of 
being  done,  and  prevents  the  temperature  of  the  gas  from  decreasing  as  much  as  it 
would  ideally,  i.e.,  if  no  friction  were  present. 

Entropy  is  defined  in  terms  of  reversible  processes.  If  a  system 
Involved  in  a  process  undergoes  a  change  from  state  PQ,  Tc,  PQ  to  Ff,  Tf,  ?{, 
the  entropy  change  of  the  system  can  be  calculated  by  choosing  a  series  of  revers¬ 
ible  processes  whereby  the  state  variables  of  the  system  change  from  PG,  T0,  PQ 
to  Pf,  Tf,  $f.  For  the  case  of  a  perfect  gas,  this  becomes  *  dS. 

If  =0,  i.e.,  no  heat  is  added  from  surroundings,  the  process 

is  adiabatic.  If  *  0  *  the  process  is  reversible  and  adiabatic,  or 

isentroplc. 

Again,  entropy  here  is  really  entropy  per  unit  mass.  From  the  first 


T 


c  IT  _  dP 
Cp  T  R  P 


This  comes  about  by:  +  £w  «  dE,  £q  -  pdV  «  CvdT,  St',  -  d(pv)  +  vdp  =  Cydt, 

S  Q  '  d  (RT)  ♦  vdp  -  CydT,  SQ  ♦  vdp  «  Cpdt. 

Note  that  for  dv  positive,  £w  is  negative,  since  <jw  is  work  done  on  the  system. 


MAC  AttS 


liMOLASSlElEH 


y^Hot^uardi 

VAN  NVrt,  CAll'OtNtA 


■imr. 


5978 


and 


S-£ 


<4 


_p 

Po 


or 


P  * 


( 

e" 


s-s0 

IT 


i 


(15) 


IV. 


STEADY  STATE  EQUATIONS  OF  MOTION 
A.  Description  of  Equations 


In  describing  the  equations  of  motion,  ve  shall  refer  to  two  differ¬ 
ent  flow  configurations,  namely,  the  two-dimensional  case  and  the  axisymmetrlcal 
case.  In  the  two-dimensional  case,  we  use  the  Cartesian  coordinates  x  and  y.  In 
the  axlsymmetrice.l  case,  the  x-axis  is  the  axis  of  rotation  of  the  x-y  plane,  and 
flow  conditions  will  be  the  same  on  any  x-y  plane  which  has  been  rotated  at  any 
angle  with  respect  to  a  stationary  line  orthogonal  to  ehe  x-axis. 

We  shall  also  use  the  steady  state  Eulerian  formulation,  whereby 
solid  boundaries  in  the  flow  are  stationary  with  respect  to  the  observer. 

The  steady  state  equation  of  ccntlnuity  is,  for  no  sources  or  sinks. 


V  •  (?q)  =0 

Where 

f  =  Density 

q  =  Velocity  vector  =  (u,  v),  u  =  q  cos  e,  v  =  q  sin  a 
Written  out,  this  becomes 

e 


(16) 


du  .  dv 

35  +  3? 


+  u  S?  +  v  ~ 


3y  "  y 

Where 

£  =  0  for  the  two-dimensional  case 
£  =  1  for  the  axlsymmetrlc  case 

This  equation  rays  that  mass  Is  preserved  at  all  points  In  the  flow  stream. 


(17) 
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The  ga3  enthalpy  for  the  equilibrium  reservoir  gas,  for  which  q  =  ( 0, 0)  is 

y  p0 

*o  3  y  —  -  ,  or  stagnation  enthalpy.  ( 

If  i  is  constant  in  the  entire  flow  field,  the  flow  is  called  isenthalpic. 

The  equation  of  energy,  Equation  (20),  is  valid  along  the  flow 
streamlines.  This  equation  could  also  be  written  as 


i ?  +  v^ 


+  Cp  T 


Cp  T0 


(20)1 


Where  PQ  =  R  P0  T0. 

Thi6  equation  expresses  the  interchangeability  of  random  and  uniform  energy  of 
motion.  The  term  C  ' T  is  related  to  random  motion,  and  the  term  l/2  (u^  +  v^)  is 
related  to  uniform  motion.  If  the  flow  is  isenthalpic,  the  energy  equation. 
Equation  (20),  is  valid  at  all  points  in  the  flow  field,  and  not  Just  along  a 
particular  streamline,  since  iQ  is  constant. 

The  temperature  T  is  that  which  would  be  measured  by  an  observer 
moving  along  the  streamline  with  speed  +  v 3?  +  v^  . 

The  perfect  gas  law  ,  Equations  (l)  or  (2),  also  holds  for  quSCo 
with  a  net  velocity.  However,  in  applying  it,  we  must  understand  that  the  quan¬ 
tities  P,  T,  and  P  are  those  which  would  be  measured  by  an  observer  moving  along 
with  the  gas  at  the  same  net  velocity,  i.e.,  they  are  static  quantities. 

Define  the  quantity  a  to  be 


We  will  later  show  that  this  is  the  local  speed  of  sound  in  a  gas.  For  a  perfect 
gas,  this  becomes 


(22)1 


since  for  constant  entropy, 


PaC^’4»  =  cyP 


HPOfit, 
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Another  relation  which  we  will  need  is  Crocco's  result, 
q  x  (\/  x  q)  =  GRAD  i  -  T  GRAD  S 


£.  Special  Results 


We  shall  write  down  seme  relations  which  will  be  needed  later. 
From  Equations  (22),  (23),  and  (3,5), 

dP  p  a2  dS  3 

^  Si  *  &  Si 


(^a) 


dP  pa2  dS  2  df5 
Si  =  pr  3y  a  5y 


Where 

>1  a 

s  = 

cv 

C.  Results  for  Isentrcpie,  Isenthalpic  Flow 

From  Equations  (22),  (22)1,  (20),  and  (2l),  and  by  differentiating 
Equation  (20),  and  by  using  Equations  (2ta)  and  (20b)  with  the  assumption  that 
S  =  constant,  we  obtain 

„  du  dv  a2  df*  _ 

+  <25a 


u  ►  V$1  *  a2  dP  a  0 

U  r  V*r—  -f  — X""  a 

oy  dy  f  dy 


i0  =  Constant  In  the  entire  flow  field. 


By  substituting  Equations  (25a)  and  (25b)  into  the  continuity 
equation.  Equation  (17),  we  obtain 


[$*■£]-$  [“s 


f  uv  ♦  UV  & 

Ox  oy  ox  oy 


MfAAt 


397b 


Differentiating  with  reapect  to  x 


du  ^  5v  .  a  f  y  p  ] 

+  v  Si  Si  [7^1  Fj  "  0 

y  ft]  p  y  f  1  3p  P  del- 1  f  y  Jap  p  apj 
5S  [  irrTj  T  ’  r?T  [T  ^  ^  3SJ  -  t  "  F  ssj 

> 

■  i  <xh>  [■’  -  f]  M  -t<rh> 

■  7 [?^ij  11 ' t  "7  <M> 


We  can  now  write  down  the  final  equations  for  isentropic,  isenthalp- 
io  ccr.cressible  fluid  flow  with  unknowns  ^,and  ’ 


From  Equation  (2.3), 


dv  du 
3x  ”  3y 


0 


(27) 


Pevritir.g  Equation  (2 6), 


du 


dv. 

3?J 


4—  +  uv  4^ +  uv  4~ 

a2  ox  oy  dx 


■f 


(26) 


Also  ve  can  write 


§;  dx  +  ^dy  -  du  (28) 

4^  dx  +  4^  dy  -  dv  (29) 
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D.  Results  for  Nonlsentroplc  Flow  (Rotational) 

Substitut.  ig  for  and  from  Equations  (24a)  and  (24b)  in  the 
ox  oy 

Navler-Stokes  Equations  (l8)  and  (19)<  we  obtain 


du  du  a2  dS  .  a2  dC  n 

u^  +  v^  +  r^  r^  “  0 


u  dv  ..  dv  +  a2  d§  +  a2  df  Q 
U  chT*  '5y  ° 

Rewriting  the  equation  of  continuity. 


Also, 


du 


dv 


de 


dS 


(30) | 

(31) 


+  u  II *  v  =  * 6  (^) 

(17) 

^  dy  -  du 

(32) 

dv  .  . 

3y  dy  "  dv 

(33) 

-  dt 

(34) 

<3 

II 

« 

(35) 

The  condition  that  entropy  be  constant  along  a  streamline  except  for  passage 
through  a  strong  shock  is 


ds  ds 
+  v3? 


-  o 
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Note  that  the  equation  of  energy  is  net  included  ir.  this  set  of 
equations.  The  following  case  of  flow  with  rotation  is  considered  here  --  there 
are  variations  in  the  rest,  enthalpy  (  i0)  and  in  the  entropy  (a),  from  streamline 
to  streamline.  The  entropy  is  a  constant  alcng  a  particular  streamline  before 
and  after  the  shock,  while  there  is  a  discontinuous  increase  in  entropy  on  each 
streamline  in  crossing  the  shock  front. 

E.  Weak  and  Strong  Waves 

The  plug  nozzle  program  (Section  V  below)  carries  out  the  solution 
cf  Equatlens  (27),  (26),  (28),  ar*.1  (29)  by  a  numerical  integration  using  the 
method  of  characteristics.  Before  going  into  this,  however,  it  is  advisable  to 
consider  the  physical  basis  for  the  mathematical  development .  Indeed,  this  is 
really  necessary,  because  certain  initial  conditions  a r.d  boundaries  .oust  be 
prescribed  before  the  solution  of  the  equations  can  be  carried  cut,  and  these 
are  obtained  from  physical  cobs iderat ions . 

We  begin  by  writing  the  continuity,  Navier-Stckes,  and  energy 
equations  for  one  dimension,  i.e.,  the  x-direction,  in  differential  form. 


du  ,  dP 

T  T  " 

0 

(36) 

udu  0 

0 

(37) 

udu  +  C  dT 

P 

-  0 

(38) 

sketch : 


F.  Speed  cf  a  Normal  Weak  Wave 

Assume  that  we  have  a  flow  configuration  as  shewn  in  the  following 


P,  P 


u 


p  +  Sp,  P  +SP 

- ^  u  +5u 


Assume  that  there  exists  a  disturbance  normal  to  the  flew  velocity  u  such  that 
the  variables  P,  P  ,  and  u  undergo  small  changes  S  P,  S  P  >  and  Su,  where 

<  <  i,  £»  <  <  1  and  such  that  the  disturbance  is  stationary. 


S_P 

P 


<  < 


If  these  changes  are  sufficiently  small,  we  may  substitute  SP  for 
aP,  SP  for  dp  ,  and  $  u  for  du  in  Equations  (.56),  (57),  and  (58)  with  a  small 
negligible  error  resulting. 


uiTOOTFmr 


- 15 
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Combining  Equations  (37)  and  (38) 


dP  „  *  .  ,  P. 

F  =  cpdT  =  T7!  d  (F> 

so  the  change  in  entropy  <s  negligible.  Combining  Equations  (36)  and  (37) 


dP 

d? 


(39* 


(Uo) 


Therefore,  for  such  a  disturbance  to  exist,  the  flow  velocity  u  must  be  given  by 
Equation  (40). 

If  we  Impose  a  velocity  u  to  the  left  in  the  above  sketch,  the  dis¬ 
continuity  advances  into  undisturbed  fluid  with  speed  u.  The  most  common  ex¬ 
ample  of  such  a  wave  is  a  sound  wave,  and 


is  called  the  local  speed  of  sound. 


C.  flanklne-Hugcnlot  Equations  for  Normal  Strong  Waves 


sketch 


Assume  that  we  have  the  flow  configuration  shown  in  the  following 

->*2 

“2  -  U1 


Pl'  / 


~>ul 


P2,  ?2 


Assume  that  there  exists  a  stationary  disturbance  such  that 


U1 


Pa  -  pi 

pi  ' 


-?i 


are  not  small  compered  to  one. 


We  refer  to  this  disturbance  as  a  strong  wave  or  shock  wave.  The 
Ranklne-Hugcniot  shock  conditions  are  now  expressed  as  Equations  (til),  (^P),  and 
(43).  The  mass  flow  across  the  wave  must  be  the  same  as  behind  the  wuve,  so 


Ug  *  a 


(J‘D 
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The  increase  In  momentum  cf  the  gas  per  unit  tine  must  equal  the 
r.et  force  on  the  gas  in  the  same  direction,  so 


-  ?2  +  ?1 


(u2  -  ux)  =  U2~  '  ui 


By  conservation  of  total  energy  (uniform  and  random),  and  using  the  energy 
equation, 

* -X_ !i  =  +  1_l i  (a*)2  (^) 

2  y  - 1  er  2  M(“2 

where  a*  -  the  speed  wherein  u  and^a  are  equal. 

The  basic  equation  for  the  velocity  change  across  a  normal  shock 
wave  may  be  derived  as  follows.  Combine  the  energy  equation,  Equation  (bj),  and 
the  momentum  equation,  Equation  (;-2),  to  obtain 

,  J  *  +  1  a*  .  V  -  1  1  ,,  ,  . 

ui '  °2  “  (ui  -  ^\~rr  ^  ~  (u) 


The  solution  to  this  equation  for  ^  is 


U1  °S  *  a* 

We  can  conclude  from  Equation  (^5)  that  if  uj_  is  greater  than  the 
speed  of  scund,  then  u2  must  be  less  than  the  speed  of  sound. 

Define  —  =  M  =  The  Mach  number  for  compressible  flew. 

H.  Mach  Waves 

Assume  that  we  have  a  weak  wave  inclined  to  the  direction  of  flow 
behind  the  wave,  which  is  stationary,  as  shewn  in  the  following  sketch. 
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Assume  also  that  di^  Is  very  small  compared  with  u^.  This  is  defined  to  be  a 

Mach  wave.  We  know  from  the  study  cf  a  weak  normal  wave  that  the  speed  of  flow 
normal  to  the  wave  must  be  a.  From  the  flow  diagram, 

•*»?-  ^  -  £  <“s> 


From  Equation  (.-6),  ve  may  conclude  that  uj_  must  be  supersonic,  and 
that  there  Is  only  one  angle  ^  at  which  the  wave  may  be  inclined  for  initial 
speed  uq. 

With  reference  to  the  bent  wall  boundary.  If  tie  flo!  Is  required 
to  he  parallel  to  the-  wall,  then  the  bend  creates  the  disturbance  which  produces 
the  wave. 

The  particular  wave  shewn  here  Is  an  expansion  wave.  If  the  bend 
had  been  upward,  It  would  have  been  a  compression  wave,  i.e.,  the  pressure  In¬ 
creases.  This  can  be  seen  from  the  relations 

d8  .  1  Jdun  cos/j, 

and 

de  sin  ^  s/m2  -i  ~ 


$  u  du,j  sin/? 


From  the  first  relation,  If  d3  Is  as  shown,  du  >  0,  from  the  second  relation, 
dP  <  0. 

I.  One-Dimensional  Nozzle  Flow 

We  new  derive  a  result  for  Mach  number  M  where  variable  cross- 
sectional  area  la  present. 


dP 

P 
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First  vc  write  the  continuity  equation  to  take  into  account  the  area  variation. 

(1*7) 


dP  du  dA 

+  —  +  —  =0 
p  u  A 


Rewriting  the  Kavier-otukes  equation. 


udu  +  a  0 


T 


Equation  (37)  may  be  rewritten  as 


udu  +  a2  =  0 


Elir.inat  ing  dP/p5  between  Equations  (37)^  and  (U7), 


du  (1  .  M2)  =  .  dA 

u  A 


(37) 


(37)1 


(*»8) 


If  dA/A  =  0,  there  are  two  possibilities:  either  M  =  1  or  du  =  0.  If,  however,] 

M  =  1  somewhere  In  the  flow,  then  at  that  point  dA  =  0.  C uch  a  point  i6  calleo 
a  throat  (local  minimum  in  the  crcss-sectional  ai-ea  A). 

V.'e  now  investigate  the  muss  flow  through  the  ccnf iguratlon  shown  in  j 
the  above  sketch.  Take  P  to  be  the  external  pressure. 


f  u  A 


2 

p  T 


V-  i 


S-  O _  p  f  x  \°  , 

-  1  Po  (pj  1  ■ 


The  mass  flew  m  has  a  maximum  at  the  point  dm/d  ~  =  0,  or  at 


VC 


-L  -  (  it  -) 

PQ  *  +  1 

For  channel  flow,  we  write  dewn  results,  assuming  entropy  constunt. 
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h  a.  X-_ 

?  v  ^ 


>f) 


rht 


(51.) 


Equation  (51)  cones  directly  from  the  energy  equation.  If  ve  put  M  »  1  Into 
Equation  (5 l),  we  get  Equation  (JO).  The  channel  outlet  is  the  point  dA/A  =»  0. 
Therefore  at  P  correspond  l  n.q  to  dm/d  —  -»  0,  M  ~  1  at  the  outlet. 

“o 

Now  assume  that  P,  the  external  pressure,  is  reduced  below  this 
value,  The  mass  flow  cannot  decrease,  as  equation  (^9)  would  predict,  when  ve 
lever  P,  but  neither  can  it  Increase.  We  still  have  M  =  1  at  the  outlet,  and 
?cut let  given  by  Equation  (50).  Therefore  we  ore  led  to  the  conclusion  that 

P  /  Poutlet  for  p  <  P  given  by  Equation  (50) 

This  means  that  there  is  a  pressure  discontinuity  at  the  outlet  (or 
shock  wave)  In  this  case.  We  may,  hewover,  add  an  expanding  section  to  the  nozrlej 
to  reach  lower  pressures  beyond  the  throat.  We  will  then  have  a  convergont- 
uivergent  nozzle. 

J .  Ranklre-Hugor.lot  Equations  for  Oblique  Strong  Waves 
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Contiuuity i 

Cp  Uj_  sin  @  ~  (ag  sir  (3  -  v2  cos  $  ) 


(52) 


Conservation  of  momentum  normal  to  the  shock  wave: 
Pj  +  P^  sin*"’  p  =  Pg  +  Pg  (a,  sin  P  •*  Vg 


Conservation  of  momentum  parallel  to  the  shock  wave: 

P 1  uj_p  sin  P  ccs  =  ^2  (u2  sin  $  *  v'2  ocs(?)  (ug  ccs^+  vpsin  $  ){5!*)i 


Conservation  of  energy  across  the  wave: 

1  ?.  .  y  P1  (J  +1)  ,2  1  ,  3  2V  '4 

C  U1  +  7^~  IT  =  a*  "  2  (u*  +  v2  >  +  *“ 


P2 

£ 


(55) 


wave. 


"he  enthalpy  will  change  when  the  flew  encounters  an  oblique  strong 


We  Just  mention  the  existence  of  reflected  waves,,  which  will  occur 
if  a  wave  is  oblique  to  a  straight  boundary.  The  reflected  wave  keeps  the  flow* 
parallel  to  the  vail. 

In  the  preceding  development,  we  have  classified  disturbance  waves 
which  can  exist  in  comprrssihie  fluids.  For  our  immeuiate  purposes,  the  class 
of  weaic  va/cs  is  most  important.  We  are  now  in  a  position  to  interpret  the  re¬ 
sults  of  tne  numerical  integration  to  follow. 

K .  The  Method  of  Characteristics  for  Isentroplc,  Isenthalplc  Flew 

By  neglecting  viscous  effects  (both  in  the  stream  end  at  the 
boundaries),  small  changes  in  entropy  and  enthalpy  across  weak  waves,  and  essen¬ 
tially  assuming  that  Cv  and  Cp  are  constant,  we  arrived  at  Equations  (27),  (26), 
(28),  and  (29).  further  assume  that  no  strong  waves  will  exist  in  the  flew 
field.  If  they  should  occur,  their  position  and  shape  would  have  to  be  deter¬ 
mined  and  the  Runklne-Hugonlot  equations  used.  In  other  werds,  strong  waves  must 
in  general  be  specified  as  "boundary  conditions"  in  tne  flow  field. 
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|  Rewriting  Eq  ations  (27),  (26),  (28),  and  (29)/ 

(27) 

(26) 

dx  +  dy  =  du  ( 28 ) 


^u  5v  .. 

3y  ■  ’Si  3  0 


cH\  dv 

? 

^  +Ty 

a2 

X* 

uc  -f 


3x 


5u  dv  2  dv 

S  +  uv  Si  +  v 


5y 
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tj  dx 


^<i,v  - 


dv 


(29) 


We  shall  refer  to  these  equations  shortly. 


Making  the  assumptions  which  we  have  mentioned,  we  note  that  «e  are 
dealing  with  systems  of  partial  differential  equations  which  are  first  order, 
that  is,  the  highest  order  derivatives  which  appear  are  first  derivatives.  (If 
we  Included  viscous  effects,  we  would  have  a  second  order  system  or"  equations.) 
This  means  that  vn  have  to  solve  an  Initial  value  problem,  l.e..  we  specify  1 

boundary  conditions  on  only  part  cf  the  boundary  of  the  flow  field.  We  cannot 
in  ceneral.  have  a  closed  boundary  px’oblem. 

In  deciding  which  numerical  scheme  to  use,  there  exists  the  possi¬ 
bility  of  writing  difference  approximations  to  our  system  of  equations  using  a 
rectangular  net.  In  order  for  these  difference  approximations  to  be  valid,  we 
require  continuity  of  higher  order  derivatives  so  that  the  mean  value  theorem  can 
be  applied  tc  the  Taylor's  series  to  obtain  a  remainder  term  which  expresses  the 
size  cf  the  truncation  error.  For  accurate  solutions,  we  of  course  require  that 
the  truncation  error  be  at  least  or  order  of  magnitude  smaller  than  the  solutions 
themselves. 


In  the  case  cf  plug  nczzle  flow  with  uniform  initial  conditions, 
however,  it  is  evident  that  if  we  assume  continuity  of  the  first  derivatives 

\  \  V  V 

and  ^  we  can  obtain  only  uniform  flow  throughout  the  flow  stream  by 
the  usual  processes  of  numerical  integration  using  rectangular  nets. 
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Consider  the  case  of  nozzle  flew  with  uni f era  initial  conditions 
prescribed  along  a  non-characteristic  line  in  a  region  cf  supersonic  flow.  Equa- 
tiens  (26)  and  (27)  are  then  a  set  of  hyperbolic  equations 


Initia 

Line 


By  the  Initial  value  theorem,  the  flow  in  a  three-sided  region  in¬ 
cluding  the  initial  line  will  also  be  uniform  (See  diagram  above).  The  lines  1 
and  2  are  characteristic  lines.  The  flew  outs’ide  of  this  three-sided  region  is 
nenuniferm,  because  of  the  curvature  cf  the  boundary,  and  the  fact  that  the 
boundary  13  a  streamline  of  the  flew.  Thus  the  flow  changes  from  uniform  to  non- 
uniform  across  the  characteristic  lines  1  and  2.  Since,  as  will  be  indicated, 
the  derivatives  at  a  characteristic  are  discontinuous,  a  valid  Taylor’s  series 
with  small  remainder  cannot  be  present  at  the  characteristics.  A  method  must  be 
used  whereby  solutions  cf  the  differential  system  cf  equations  can  be  patched 
together  along  characteristic  lines. 

Courant  pcir.t3  out  that  the  nonlinear  equations  cf  hydrodynamics 
(i.e..  Equation  (27))  exhibit  this  behavior  of  discontinuities  being  generated  in 
the  flow  field,  even  though  ncr.e  are  present  cn  the  initial  line.  This  cannot 
happen,  however,  for  linear  sets  of  equations  of  hyperbolic  type;  for  example, 
the  wave  equation 


2  ,o2u 

* 


d2u 

+ 


Ve  Just  mention  that  one  vay  of  eliminating  the  characteristics  is 
to  a  second  order  viscosity  U  rm  to  the  energy  and  momentum  equations. 
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This  second  order  term  has  little  effect  except  at  places  where  the 
derivatives  ux,  Uy,  vx>  Vy  otherwise  become  discontinuous  (at  characteristics), 
then  it  has  the  effect  of  smoothing  out  the  transition  lines  so  that  the  charac¬ 
teristics  vanish. 

One  method  of  finding  the  characteristic  lines  for  Equations  (27), 
(26),  (28),  and  (29)  is  now  illustrated.  We  can  write  these  equations  in  the 
following  form: 


pd  -4) 

ac 

-  4  uv 

-  4  uv 

(1-4) 

C h} 
dx 

- 

y 

0 

1 

-1 

0 

du 

0 

dx 

dy 

0 

0 

dv 

35 

du 

0 

0 

dx 

dy 

dv 

dy 

• 

dv 

(56) 


Solving  for  —  by  Cramer's  rule,  we  obtain 
dx 


du  (a2  -  v2)  dv  dy  -  2  uv  dy  du  -  (a2  -  v2)  dx  du  +  €  a2  -w-  dy2 

"5x  2  2  2»  , "  ("2  ( 57 ) 

or  (dx^  +  dy*)  -  (udy  -  vdx) 


du 


For  • —  to  be  indeterminate,  it  is  necessary  and  sufficient  that  both  numerator 
dx  ^ 

and  denominator  be  zero.  This  implies  that  —  ,  r—  and  —  are  also  indetermin¬ 
ate.  dy  dx  dy 

Discussion:  Since  the  denominator  is  zero,  the  coefficient  matrix  is 

singular,  hence  of  rank  les6  than  four.  This  implies  a  linear  relation  among 

,  ^,and  ar.d  hence  and  are  indeterminate, 

dx  Jy  dx  dy  oy  dx  dy 
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Verification:  Assume  Equation  (62)  Is  true.  Then 


uv  +  a  Vq^  -  (0  +  c^)  -  *an  2.  1  tan  tan  9  +  Vq^  -  a^ 


2  2 
u*  -  a 


1  -  tan  9  tan<  1  -  tan  Q 


o2  sin  3  ccs  3  +  b2Vm2  -  1  ?  Vm2"  -  1  sin  9  ±  cob 


2  -2  2 

q  cos  **  - 


Vrf 


1  cos  9  -  sin  9 


q2  sin  9  cos2  oVm2  -  1  +  a2  (M2  -  1)  cos  9  -  .q2  sin2  3  cos  3  -  a2  sin  s'/m2 


=  q2  cos2  3  sin  9Vm2  -  1  +  <q2  cos^  9  -  a2VM2  -  1  sin  9  -  a2  cos 


a2  (M2  -  l)  cos  9  -  q2  sin2  3  cos  3  =  q2  cos^  9  -  a2  cos  9 


q2  -  q2  sin2  9  £  q2  cos2  9 


2  2 
9  -  q  =  ° 


Similarly,  we  could  verify  Equation  (63). 

The  numerator  of  the  quotient  of  Equation  (57)  is  set  equal  to  zero.  It  is 


(a 2  -  v2)  dvdy  -  2  uv  dydu  -  (a2  -  v2)  dxdu  +  fr  a2—  dy*-  =0  (6 

y 


The  denominator,  set  equal  to  zero.  Is 


(a2  - 


(6 


u2)  dy2  +  2  uv  dxdy  +  (a2  -  v2)  dx2 


0 
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Multiplying  Equation  (6U)  by  dx  and  substituting  (a2  -  v2)dx2  from  Equation  (65), 
we  obtain 


(a2-v2)  dvdydx  ■  2  uv  dydudx  +  (a2-u2)  dy2du  +  2  uv  dxdydu  +  fa2  ^  dy2  dx  *  0 


(a2-v2)  dv  dydx  +  (a2-u2)  dy2  du  +  6  a2  ^  dy2  dx  »  0. 

Ignoring  the  possible  root  dy  =  0,  we  finally  obtain 


(a2-v2)  dvdx  +  (a2-u2)  dy  du  +  £a2^  dy  dx  ■  0 


Substituting  Equation  (58)  in  Equation  (67),  and  ignoring  the  possible  root  dx 
0,  we  obtain 


(a2-v2)  dv  +  (a2-u2)AL  du  +  £a2  ~  dy  *  0  (66) 

Substituting  Equation  (59)  in  Equation  (67),  and  Ignoring  the  possible  root  dx  » 
0,  ve  obtain 


(a2-v2)  dv  +  (a2-u2)AR  du  £a2  -  dy  »  0 

y 


Equations  (68)  and  (69),  give  relations  between  u,  v,  and  y  which  must  hold 
along  the  characteristic  lines  given  by  Equations  (58)  and  (59).  These  are 
referred  to  as  compatibility  equations. 

We  will  now  go  through  a  developinent  to  put  Equation  (68)  in  a  mere 
convenient  form. 
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From  Equation  (.16), 


dT  = 


Frcra  Equations  (77)  and  (75), 


i— [ 

M  ..2, 2  l 


(1  +  ±L±  m2) 
2 


(  K-l)  M 


dM 


= 


Tp  Cp 


^  (i ,  ta  „2, 


(K-l)MdM 


Since  q  •=  Ida  and  a  =  s/ JT  RT, 


To  cp 


la 

a 


(i  *  m2) 

2 


-5  (  /  -1)  M  d  M 


2T  RM 


£~  (X  -  l)dM 


2  y 

M  i  R 


1  + 


t-  1 


M2 


(1  +  M2  ) 


dM 


M  ( 1  +  1^-1  M2) 
2 


From  these  results,  Equation  (71)  assumes  the  form 


d«  _+C| 

M  ( 1  +  ^  7  1  M2 ) 


2 


tan  8 


V'm^t 


1  -  ton  9 


]dr. 


597a 


=  0 


(corresponding  to  X  L) 


(77) 


(78) 


(79) 


(7l)J 
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The  unique  solution  to  Equation  (82)1  i8 

9  =  -<*,  (Left  line),  (84) 

and  the  unique  solution  to  Equation  (83)^  Is 

f 

9  =  .  (Right  line).  (65) 

To  Interpret  equations  (84)  and  (85),  let  us  refer  to  the  following 

diagram: 


If  9  =  -  °<  for  the  left  line,  then  obviously  the  left  line  lies  along  the  x-axis, 

If  9  =  for  the  right  line,  then  the  right  line  lies  along  the  xaxts. 

Me  might  comment  here  that  the  notation  "left  characteristic  line", 
or  "left  line”  denotes  the  direction  that  a  ball  would  roll  If  placed  along  the 
line.  Likewise,  a  ball  would  roll  to  the  right  if  placed  on  a  'right  line". 


As  a  practical  matter,  we  can  avoid  the  condition  dx  -  0  by  requir¬ 
ing  M>1  in  the  calculations.  This  point  will  b«..  clarified  later.  The  condition 
dy  =  0  usually  does  not  occur,  it  if  does,  there  is  s  special  provision  which 
can  be  carried  out  in  the  calculations. 

The  condition  y  =  C  Is  avoided  in  the  calculations  as  vill  be  seen 

later. 

V.  PLUG  NOZZLE  PROGRAM 

A .  General  Relations 

In  the  plug  nozzle  program,  we  deal  with  flow  which  Is  Is*  tropic 
cnJ  Isenthslpie.  This  implies  that  in  the  flow  field  there  are  r.c  vort.,city, 
heat  transfer,  or  viscosity  effects,  and  no  chocks.  The  chornoterictK which 
apply  to  this  type  cf  flow  are  expressed  by  Equations  ,  0?)  and  (S'): 
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dy  -  t«A  (6  +  °()  dx  ■»  0  (Left  line)  (62) 

dy  -  "tan  (0  -  °<)  dx  =  0  (Right  line)  (63) 

The  compatibility  relations  which  must  hold  along  the  characteris¬ 
tics  are 

d9 


tan  9 


V m2  - 


tan 


dx 

---- 

•> 


*  0  (Along  left  line)  (n)1 


d9  + 


Vm^- 1  dH 

M  (1  -r  — —  >£) 
2  > 


tan  0 _ 

-  1  +  tan 


0  (Along  right  line) (72) 1 


We  now  make  a  slight  modification  to  Equations  (71)1  and  (72)^.  In  the  existing 
program,  we  add  a  snail  correction  factor  to  the  coefficient  of  dM  in  Equations 
(71)1  and  (72)  ,  so  that 


-  -J. _ 

+  M  (1  +  ^  M2) 


is  replaced  by 


-  ^  at  “  1 _ 

+  M  ( 1  +  M^) 


1 


"<  ('.«)  •  57  * '  m  I 


This  comes  about  Ly  assuming  that  Cp  and  Cv  are  not  constant,  but  ^  »)C(T). 
Wt  otlll  use  the  perfect  gas  law  as  our  equutlon  of  state,  however. 
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Derivation  of  the  new  expression  for  the  coefficient,  of  dM- 

Rcc.  11  that  the  dM  term  was  originally  i  ^  (  Equations  (71) 

and  (72)).  q 

By  the  equation  of  energy  in  differential  form, 

CpdT  +  q  dq  =  0  (Along  a  streamline).  (88) 

This  can  be  rewritten  as 


We  have 


(R_L)  ±t  _  _  d* 

-  V  q2  q 


^  RT 


M2  = 


7  RT 


Equation  (88)1  becomes 


dT  = 


X  RM2 


2  q  dq  _  2q2dM  _  q2  d 't 
K  RM2  /  RM?  X  2RM2 


la 

a 


±  (RJL) 

a2  l2f- 


f-LlJs 

ms 


2qg  dM  .  qg  d/ 


^  RM5  2  RM2 


11 


(88)1 


(88)H 


Froir  Equation  (88) 


we  finally  obtain 


<<  X 
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In  line  with  this  assumption  that  Z-:  -  CV(T)  and  Cp  -  C^(T),  Equa¬ 
tion  ( 51 )  loses  significance,  since  it  is  derived  assuming  Cv  and  Cp*a.re  con¬ 
stant.  An  expression  for  dT/cW  which  can  be  derived  from  the  energy  equation  is 


(1  +  Y*  '  (T))  + 


However,  we  still  use  Equation  (51)  to  obtain  initial  line  temperatures.  Also, 
we  neglect  certain  other  effects  which  should  modify  cur  initial  assumptions. 

The  general  type  of  nozzle  which  will  be  used  is  illustrated  in  the 
following  diagram,  which  represents  a  cress  section.  In  the  axisymmetrical  case, 
this  cross  section  is  any  plane  which  has  included  the  x-axls,  or  the  axis  of 
symmetry.  In  the  two-dimensional  case,  the  cross  section  is  any  plane  parallel 
to  the  x-y  plane.  The  lower  half  of  the  nozzle  Is  hot  shown,  since  it  is  a  mir¬ 
ror  image  of  the  upper  half. 


Outer  Wall 


y-axis 


FLOW 

FIELD 


Inner  Wall 
(Plug) 


=  £  between  outer  wall  and  x-axis  at 
the  initial  line 

*  £  between  plug  wall  and  x-axls  at 
the  initial  line 

*  Initial  guess  for  final  streamline 
angle 


axis  of  symmetry  (x-axis) 


The  variables  x,  y,  M,  and  9  are  presecrlbed  at  discrete  points  on 
the  initial  line  as  part  of  the  starting,  cr  input  data.  The  initial  line,  which 
need  not  be  vertical,  is  generally  assumed  to  be  in  the  throat,  or  point  of  mini¬ 
mum  cross-sectional  area,  and  so  M  along  the  initial  line  is  usually  assigned  a 
constant  value  Just  larger  than  one.  A  table  cf  (x,  y)  values  describing  the 
cuter  wall  and  the  plug  wall,  al6c  the  angles  Scuc  and  9ln  are  also  needed. 
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It  is  assumed  that  a  continuously  differential  curve  can  be  passed 
through  the  outer  and  inner  wall  points,  starting  at  x  =  0.  (From  the  diagram, 
it  is  evident  that  this  is  not  possible  at  point  c.)  We  will  require  also  that 
the  streamlines  of  flow  are  parallel  to  the  inner  and  outer  wall  at  points  where 
the  derivatives  are  continuous.  At  point  c,  we  must  handle  the  situation  differ¬ 
ently.  The  methods  of  procedure  for  insuring  parallel  flow  at  the  boundaries  and 
for  handling  conditions  at  point  c  are  now  described, 

B.  Physical  Considerations  at  Outer  and  Inner  Wall  Boundaries 
1.  Reflected  Mach  Waves 

We  can  provide  a  physical  basis  for  our  boundary  conditions  by 
showing  that  the  characteristics  line  segments  are  also  the  Mach  wave  lines. 
Equation  (16)  describes  the  angle  of  inclination  of  a  Mach  wave  to  the  flow 
direction  as  a  unique  function  of  M. 


A* 


sln  P  *  ir 


For  the  characteristic  line,  we  have 

.  8  1 

sin  c(  =  —  =  tt-  , 

q  M  ’ 

where  <X  =  Angle  between  the  characteristic  line  und  the  streamline. 


Obviously, 


cX  =  P 


Note  that  we  are  dealing  with  stationary  lines  here. 

Therefore,  the  characteristic  lines  are  also  Mnch  Hi  s,  i.e., 
the  loci  of  Mach  waves. 


lC 


From  physical  considerations,  when  a  Mach  wave  is  incident  on  a 
smooth  boundary  oblique  tc  the  direction  of  flew,  there  must  be  o  reflected  wave 
vh'th  renders  the  direction  of  flow  again  parallel  to  the  boundary,  as  shown 
be lew. 


Initial  Flew 


Incident 

Wave 


Final  Flow 


Reflected 


In  the  plug  nozzle  program,  the  reflected  wave  will  be  manifest  ss 
a  characteristic  line  originating  at  the  boundary  point  where  the  incident  char¬ 
acteristic  line  intersects  the  boundary.  The  point  of  intersection  is  determined 
by  the  calculations. 

2 .  Prandt 1 -Meyer  Expansion 

In  the  plug  nozzle  program,  characteristic  lines  are  generated 
starting  at  the  initial  line  and  moving  to  the  right.  Let  us  consider  what 
happens  when  we  have  reached  point  c  on  the  plug. 


Outer  Wall 


x-axls 


is  the  flrot  left  line  tc  reach  point  c. 


By  first  considering  that  theie  is  u  ir.e  curvature  at  point  <; 
and  then  taxing  the  limiting  case  of  a  corner  at  point  c,  it  becomes  evident 
that  there  must  be  a  fan  cf  characteristics  originating  at  point  c  and  cer.tino 
:.o  until  the  final  streamline  angle  3g  is  reached.  Pictorlally, 


-  it  - 


I 
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The  expansion  angle  Is  denoted  by  f?  . 

How  consider  what  happens  at  the  right  characteristics  in  the 
fan  as  we  approach  point  c.  Evidently  dx  approaches  zero.  Equation  (72)1  then 
becomes 


,  ^:7L« -  [l  +  M0t(T.  M)1 

M  (1  +  M2)  L  ■* 


and  the  right  characteristics  converge  to  the  point  c. 

Equation  (86)  gives  us  a  relation  between  streamline  angle  and  Mac  number  in 
the  vicinity  of  point  c. 


Since 


We  can  obtain  a  relation  between  9  +  °<  and  M  also. 


c<  =  sin'1  -jjj- 


mVm2”  - 


d(0  t«()  »d8  +  dd=  -  (i  +Mo<.(T,M)  +  ~~~z=  (87) 

M  (1  +  ^3^2)  M  VM2-  1_ 


M  (1  ♦  *^2) 


From  Equation  (87);  we  can  calculate  M  for  the  characteristic  left  line  ;rlgln- 
ating  at  c  at  angle  9  +  o<  by  numerical  Integration  of  Equation  (>'’■*  the 

program,  we  specify 

d(i  «-o<)  -  A(0  +  o()  =  A  ^ 

to  be  a  small  negative  constant  angle. 
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Note  that  since  3  becomes  negative,  the  "ball  rolling"  definition  of  left  line 
becomes  invalid!  Suffice  it  tc  say  that  the  left  lines  which  are  dealt  with 
here  retain  the  property  that  the  angle  is  added  to  the  angle  3  to  obtain 
the  angle  of  inclination. 

C.  General  Description  of  the  Catculative  Procedure 

By  referring  to  the  program  abstract,  we  can  determine  what  input 
data  are  required  for  a  calculation. 

1.  Curve  Fitting 

The  first  procedure  is  to  determine  continuously  differentiable 
functions  which  pass  through  the  specified  cuter  and  inner  wall  points.  The 
slope  of,  say,  the  cuter  wall  function  at  both  ends  of  the  nczzle  will  be  that 
specified  by  the  input  data. 

This  curve  fit  is  obtained  in  the  following  way.  We  begin  at 
the  first  cuter  wall  point,  where  the  slope  Q0  p  is  specified.  (The  inner  wail 
calculation  is  carried  out  in  the  same  way)  Tass  a  parabola  through  the  first, 
second,  and  third  specified  outer  wall  points,  and  determine  the  slope  30  g  of 
this  parabola  at  the  second  point.  Then  determine  the  coefficients  of  a  cubic 
which  passes  through  the  second  and  first  point,  and  which  has  slope  30  p  at  the 
first  point  and  slope  ©0,2  at  the  second  point.  This  cubic  Is  the  analytic  ex¬ 
pression  for  the  cuter  wall  between  the  first  two  inp>ut  points.  Obtain  a  cubic 
for  the  second  and  third  points  in  the  way  just  described,  by  passing  a  parabola 
between  the  second,  third,  and  fourth  points,  etc. 

2.  Characteristics 

In  generating  the  characteristic  net,  we  are  primarily  ccr.cerr.ed 
with  finding  points  cf  intersection  of  left  lines  with  right  lines,  or  with 
points  of  intersection  of  left  lines  with  the  cuter  wall,  or  right  lino's  with  the 
inner  wall. 
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We  will  nov  describe  in  detail  the  method  for  calculating  T 
X3,  Y3,  M3,  and  S3  at  point  3  in  the  diagram  above,  where  we  know  Tj.,  Xi,  Y 1 
Si,  Mi  and  72,  X2,  Y2,  M2,  S2- 

Express 

dxL  =  x5  -  Xj_,  dyL  =  M3  -  yi3 

dx1*  =  X3  -  x2,  dyR  =  y3  -  y2 

Rewrite  Equations  (62),  (63),  (7l)^,  and  (72)L 

dyL  -  AL  dxL  =  0,  or  y3  -  yx  -  AL  (X3  -  xx)  -  0 

dy^  -  AR  dxK  =  O,  or  y 3  -  y2  -  AR  ( *3  -  *?)  '  0 


s3  *  ®1  '  aL  (m3  '  «l)  +  BL  (xj  -  x3)  0 


S3  -  02  +  AR  (M3  -  Mg)  -  flR  (X3  -  x2)  =  0 


i/here 


-_yZ._i.-L - Fl  +  Mo((T,  M)  j 

M  (1  +  ~~  M2)  >■  i 


AR  »  -VZH - [1  +  Md(T,  M)  1 


a*r  cm 
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g  _ _ tan  @  _ 

y  "Vfi^  ~  1  -  tan  9 


Vk2  -  3.  +  tan  9 


AL  =*  tan  (8  +  <K) 


tan  (9  ■**  o< ) 


We  also  have  equation  (5l)  which  gives  T  =  T(M). 

for 

The  star  indicates  that  the  starred  quantities  A  ,  A  ,  B  and 
a  are  to  be  evaluated  at  points  on  the  lines  connecting  Xj  with  Xg  Qn(*  *!• 
That  is,  A1^  and  are  to  be  evaluated  at 


9,  +  9. 


Mx  +  Mj 


Tj  ♦  t3 


rather  than  at,  say  9j_,  M,  and  T^. 

This  more  accurate  procedure  Is  called  the  mean  value  lattice  point  method. 

Since  this  implies  that  Equations  (62),  (63),  (71)^"  and  (72)1 
are  nonlinear,  we  cannot  solve  directly  for  X3  ,  Y3  ,  ,  and  9<,  tut  must  use 

an  Iterative  procedure. 
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.  To  begin  tiic 

T,,  onJ  AR,  AR  ,  und  IrR  at  K, 
(71)1  and  (72)^  for  the  in it  in 


To  begin  t}ic  it*1  rat  j  <  n,  evnluMte  A  ,  A  ,  and  B  at  utkI 

und  IrR  at  Ki,  -,2*  and  Tg.  lien  solve  Equations  (62),  (tj>), 
cr  the  initial  Iterates 


V,  VA  and  * A  where  1=1. 


Then  evaluate  iT'1  from  Equation  (yo)  us in/j  Am  =  *M'  -  Mq  and  Am  -  -  Mp, 

und  expressing 


Tq  +  T2  +  A'A 


/\  =  average  of  the  two  differences  obtained  from  -  Mq  and  *M^  -  M£ 

Since  Equation  (5l)  says 


calculating  T  Involves  an  iteration  also. 

i+1  5  1+1  5  1 TL  ,1  ,1  }  1  jJ.  ,  v  . 

Nov  obtain  x  and  y  from  A  =  l/<‘  tun  (  r>  +  «  ;  +  tun(d^+°lq) 


Obtain 


J^R  =  1/2  j^tun  ( -  i«x5)  +  tan  (tb,  -  oOjj 


1+V  ,  1  +  V  ,  1+V 


1  A  L  «  1/2  A1'  (V,  VA(  y  -  (T)5)  -I  A1,  (Mq,  -q,  Tq,  /  '(Tq) 


and  so  on. 


Ik'MfiW 


TF1FIT 


UNCLASSIFIED 


»i*OfT. 


« 


van  nut  a.  California 


997  d 


B 


L 


B* 


# 

tan  0 _ 

-  1  -  tan  9 


tan  9 
-  1  + 


_ 1 

tan  0  J 


* 


AL  =  tan  (8  +°<5 

A*  =  tan  (8  -o () 


We  also  have  equation  (5l)  which  gives  T  =  T(M). 

T  Q  T 

The  star  Indicates  that  the  starred  quantities  A  ,  A  ,  B  and 
s  are  to  be  evaluated  at  points  on  the  lines  connecting  X3  with  X2  and  Xi» 
Hint  is,  and  are  to  be  evaluated  at 


M  » 


3l  +  °3 


Mx  +  Mj 


and 


_  T1  +  t3 

T  = - * 

2 

rather  than  at,  say  0q,  M,  and  Tj_. 

This  more  accurate  procedure  is  called  the  mean  value  lattice  point  method. 

Since  this  implies  that  Equations  (62),  (63),  (71)^  and  ( 72 ) 1 
are  nonlinear,  we  cannot  solve  directly  for  X3  ,  Y3  ,  M3  ,  and  03,  tut  must  use 
nn  Iterative  procedure. 
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r,  :,t  l  ::u*i  this  Iteration  until 


x^-x^l  <  (f 


rv-v  <e. 


i+1M3-iM5|  <CM 


!i^l93  .  1 8  3 


For  the  case  of  outer  wall  or  inner  wall  intersection  points,  we 
;;  through  a  mod  If ication  of  this  type  of  iteration.  However,  an  cuter  wall 
Intersect : on  point  involves  only  a  left  line,  l.e..  Just  one  set  of  points  T,  X, 

'{ ,  M,  -  which  are  known.  An  inner  wall  intersection  point  involves  only  a  right 
line.  When  the  iteration  for  the  intersection  point  (x,y)  is  complete,  tan  3  will 
be  dy/dx  of  the  particular  cubic  that  represents  the  boundary  at  the  point  { x,  y). 

Generally  this  is  the  way  the  characteristic  net  is  begun. 

Assume  that  the  initial  line  is  neither  a  left  line  or  a  right  line.  Assure,  say, 
that  three  points  are  given  as  input  data  to  describe  the  initial  line. 


Cuter  Wall 


Initial 
Line  1 


Inner  Wall 


:un  calculate  T 
end  2.  (See  a' 


Ti,  X^,  Y^,  Mx,  and  9j  from  prescribed  conditions  at  Input 
above  diagram;. 
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It  should  be  noted  that  the  Intersection  points  102,  102,  or 
103j  97  do  net  as  a  rule  Intersect  the  boundaries  at  input  data  points.  Also,  we 
cannot  trace  a  streamline  by  assuming  that,  say,  a  streamline  can  be  drawn  from 
point  101,  99  to  point  102,  100.  We  must  obtain  streamlines  by  another  considera¬ 
tion.  All  we  know  in  this  connection  is  the  streamline  angle  8  at  the  two  points. 
Note  that  the  net  refinement  depends  up-cn  the  number  of  specified  initial  line 
points . 

This  process  is  continued  until  we  have  generated  a  left  line 
past  the  inner  wall  cutoff  point  c.  (To  obtain  this  line  we  must  use  an  imaginary 
extension  of  the  inner  wall  boundary  past  c.) 


We  then  Interpolate  betveen  left  lines  n  and  n+1  to  obtain  the 
first  special  left  line  which  begins  at  point  c.  Starting  with  this  line,  we 
develop  a  fan  of  left  line3  about  point  c,  evaluating  M  ori  each  line  using  Equa¬ 
tion  (87)  - 


d  (9  +<*) 


f-V^ZI -  (1  +  MX(T,M)  +  — 1 

M  (1  +  -  M2)  mVm2  -  1 


dM  (87) 


where  d  (9  +<*)  -A(9  +o0  is  a  specified  constant. 

This  expansion  is  continued  until  9  corresponds  to  -5 ,  the  creel 
fled  streamline  ungle  at  point  c.  For  this  correspondence  to  occur,  we  must  sub¬ 
divide  A(9  ♦«)  near  the  angle  9a . 

If  the  initial  guess  for  9S  is  good,  we  can  trace  a  streamline 
starting  at  point  c  which  is  asymptotic  to  the  x-axis. 

If  it  lr  net  good,  we  can  estimate  a  new  9S,  and  try  rgain. 

1 


Along  this  streamline,  Mg  and  ?g  are  constant,  therefore,  we  need  only  c c -  out 
(X,  Y,  8}s  at  discrete  pe'nts,  say  point  5  above. 
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•  In  calculating  ( X,  Y,  9)  at  point  },  we  use  the  squat  1  ;;n  to.,  the  j 
streamline  i 


dy  =  tan  9  dx 


or 


(y^  -  .V2)  3  tan  ( .  ---■?)  (*3  -  x.p) 

and  Equations  (63)  and  (72)^  for  right  line  <  1,  3  > 

Eventually,  we  will  generate  a  left  line  which  will  extend  past 
the  last  outer  wall  data  point.  When  this  happens,  we  extend  the  cuter  wall  to 
intersect  the  line,  and  then  proceed  as  Indicated  below,  we  have  omitted  the  right 
lines. 


3 .  Nozzle  Efficiency 

In  order  to  calculate  nozzle  efficiency,  we  must  first  develop 
a  system  of  equations  which  describe  an  ' ideal"  nozzle. 

To  do  this,  we  develop  briefly  a  one-dimensional  theory  which 
Ideally  might  apply  to  a  nozzle  cf  slowly  varying  crcss-noctlonal  ar-.-.i.  •  dis¬ 

cussion  follows  Hermann  and  Puckett,  Aerodynamics  of  »i  Compressible  - '  1  ul  .  In 
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From  Equations  (93)  and  (94) 


<£>'  - 


~2  [y?i  (i + Mi 


'The  stream  thrust  at  aro'  point  cf  the  nozzle  is 


F  =  pA  t  fAu2  =  pA  (ifi  X  +  l) 


Let  us  soy  that  ve  want  to  evaluate  F  at  a  particular  area  A.  If 
we  know  A *,  we  car.  calculate  M  from  Equation  (9.3)  >  once  having  determined  P  from 
Equation  (94) •  Then  we  have 


¥+  x  rr'y  2 

jf  M  *e 

y  i  ~1 

(1  +  ^  M2) 


where  y#  =  the  value  y,  the  perpendicular  distance  f rcu  the  exit  xe  to  the  wall. 
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i.’.v  we  develop  the  thrift  developed  by  the  plug  nozzle.  This 
done  at  discrete  points  x^.  corresponding  to  the  intersections  of  left 
:  12:  the  cuter  wall.  The  thrust  will  be  expressed  as  the  sum  of  four  terms 


fl  +  +  Fje  +  h 


fptie 


(  a osur.e  that  the  flow  is  leentrecie,  sc  the  stagnation  pressure  P0  is  the  same 

r.hiv  ughc.it  the  field  of  flew.  Mote  that  the  equations  apply  only  to  the  uxlsym- 
ric  case. 

The  Expression  for  Fp: 


is 


This  term  is  the  mass  flow  component  parallel  to  the  x-axls.  It 


Fi  (m2  v  + 1) 

a-. - - — — I— 


_  rr/cut,  2  in„  2 \ 

it  /o"  y0  ' 

°  u.JLJrf)** 

2 

for  flew  parallel  to  the  x-axis. 

’•■•'here  M  is  the  Mach  number  at  the  initial  point  X0  (M  is  assumed  to  be  cons  tarn, 
on  this  cross  section). 

The  Expression  for  Fg: 

This  term  is  the  force  on  the  base  cf  the  plug.  We  can  calculate 
?  in  the  static  region  from  Equation  (9-0.  since  M  is  constant  cn  the  final 
streeallr.e.  V2o  obtain 
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This  terra  exp  re  rises  the  fvreo  .  n  t  he  vuter  wall  which  odds  to 
subtracts  frera  the  thrust. 
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y  =  ye 


IX 


y  (4 ♦ 1  -  -  f  f 

Z_.  r°  *o 

.-=  0  J  V 


y  -  o 


y  =  ye 
pd  A 
y  =  C 


are*  he  unci  a  ry  points  ( outer  wall)  ana  Pj  are  preenures  at  intersections 
of  left  lines  with  the  cuter  vail.  P*  can  be  computed  from  Mj  by  Equation  (9*0 


'Hie  Expression  for  Fj( : 


This  tern  ic  sliaijiar  tc  Fj  except  it  expresses  the  force  cctr.po 
,;rit  on  the  inner  wall. 
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i  +  1 


y  =  0 


Here  Pt  are  pressures  at  intersections  of  riyht  lines  with  ihe  inner  wall.  De¬ 
fine  tmi  efficiency  to  be 


n 


^  Pile  ?o 


N  =  P0  F 


Since  we  have  ^  at  discrete  points,  v«  toiila  "bivoV  off"  the 


:uter  wall  at  one  of  those  discrete  points,  and  have  it;.*  efficiency  for  the 
;hcppo<l -off  nozzle.  How.. ver,  fnra  physical  cons iderat. ions,  we  i::us»  v.t.  chop  of 3 
:he  o viter  wall  to  the  left  if  the  right  line  which  hits  the  inner  wall  boundary 
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Symbol 

Description 

c< 

Mach  angle 

a 

Speed  of  sound 

a* 

Critical  speed 

Cv 

Specific  heat  at  constant  volume 

CP 

Specific  heat  at  constant  pressure 

E 

Internal  energy  per  unit  mass 

€ 

=  1  if  axisymmetrical  geometry,  =  0  if  two-dimensional 

* 

Ratio  of  Cp  to  Cv 

h  or  i 

Enthalpy  per  unit  mass 

k 

Boltzmann  constant 

m 

Mass  in  Section  IIT-A,  mass  flow  otherwise 

M 

Mach  number 

/ 

P 

f 

Pressure 

? 

Density,  mass  per  unit  volume 

Q 

Heat  energy  per  unit  volume 

q  Flow  velocity  =»  ( u,  v) 

q  Flow  speed  = 

R  Gas  constant 

S  Entropy  per  unit  nans 

T  Temperature 

t  Time 

9  Flow  angle 

Viscosity 

u  Component  of  velocity  in  the  x-dtrectlcn 

v  Component  of  velocity  in  the  y-dlrectlon 

V  Volume  per  unit  m ass 

W  Work 

x,  y  Coordinates  of  a  point  in  the  Euclidean  plane 
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